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Abstract 



We present results on the mixing parameter Am^ and tests of the CPT 

^~^- symmetry in the B'^B^ mixing from the time evolution of dilepton events on 

(~| ! the T{4S) resonance. The analysis is based on a 5.9 fb~^ data sample collected 

by the Belle detector at the KEKB accelerator from January to July of 2001. 

^ ■ We obtain Am^ = 0.463 ± 0.008 (stat) ± 0.016 (syst) ps"^. No evidence for 

'^ , CPT violation is found. 

a [ This is the first determination of Arud from time evolution measurements 

at T(4S'), and the first time that an experimental limit on \mno — mno\ /rrino 

d ^d d 

has been obtained. 
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1 Introduction 

The mixing parameter of the B^-B'^ system is equal to the mass difference of the 
physical (mass) eigenstates of the B^ mesons, Am^. Arrid is considered a fundamental 
parameter of the B'^ system and a crucial element in the study of CP violation. Its 
extremely small value (13 orders of magnitude smaller than the 5° mass) enhances the 
sensitivity of the B^ B^ interferometry, making the neutral meson system a powerful 
probe for physics beyond the Standard Model. 

At KEKB, 8.0 GeV electrons and 3.5 GeV positrons collide on the T{AS) reso- 
nance, producing correlated B'^-B^ pairs that live in an antisymmetric state. The 
time-dependent asymmetry between same-flavor {B^B^ /B^B^ ) and opposite-flavor 
[B^B^ ) decay pairs undergoes an oscillation with frequency equal to the mixing pa- 
rameter, Arrirf. The center of mass (CM) of the meson system is moving along the 
electron beam direction with a Lorentz boost of Pj = 0.452. Since the two B mesons 
are practically still at the T(45') frame, we can use this boost and their decay vertex 
separation, Az, to calculate the proper time difference between their decay times: 
At = Az/f3'yc. The typical decay vertex separation is 200 /im, a distance well within 
Belle's vertexing capabilities [0. 

High momentum electrons (or positrons) and muons are used for both the B flavor 
tagging and the decay vertex determination. The proper time distributions for same- 
and opposite-flavor B decay pairs can be used for the study of the 5° — i?° time 
evolution and the Arud measurement. By allowing for a possible non-zero mass or 
lifetime difference between B^ and B^ — encoded in the complex parameter cos 6 of 
the system's "effective" Hamiltonian 0] — we can use the same time distributions to 
test the conservation of the CPT symmetry in the 5° — 5^ mixing. 

The analysis presented here is based on data samples collected by the Belle de- 
tector on the T(45') resonance and 60 MeV below the peak from January to July of 
2001, corresponding to integrated luminosities of 5.9 fb~^ and 0.6 fb~^, respectively. 

2 The Belle Detector 

The Belle detector is designed to make precise measurements of charged and neutral 
particles over a wide acceptance. The combination of a silicon vertex detector (SVD) 
and a central drift chamber (CDC) surrounding the SVD is used for vertexing. For 
measurements of the charged particle momenta, the SVD and CDC are embedded in 
a 3.4 m diameter superconducting solenoid. The generated 1.5 T magnetic field runs 
parallel to the direction of the z-axis (defined by the electron beam direction). 

For the particle identification (PID), three subsystems are employed: the CDC 
(measuring the mean energy loss dE/dx), the Aerogel Cerenkov Counters or ACC (a 



set of 1188 Cerenkov radiation modules located outside the CDC volume) and the 
Time of Flight system or TOF (128 scintillation counters providing timing informa- 
tion). 

A pair of calorimeters, the Electromagnetic Calorimeter or ECL (8736 Csl/Tl 
crystals) and the Electromagnetic Forward Calorimeter or EEC (320 BGO crystals), 
is used for the detection and identification of photons and electrons. 

Finally, a set of 14 layers of resistive plate counters sandwiched between 4.7 cm 
thick iron plates and filled up with a gas mixture on high voltage is the K^ and muon 
subdetector (KLM). KLM is the only detector system placed outside the coil. The 
iron also serves as the return path for the magnetic flux. 

Ref. [0 gives details on the specifications of the detector components. 

3 Analysis 

In this analysis we exclusively use leptons for flavor tagging. This mode is referred 
in the bibliography as "dileptons": events where both B mesons are flavor-tagged 
by leptons. The term "dileptons" implies that both of the B mesons decay via a 
semileptonic process -B°(t), -Brf(i) — ^ £^X^ (where i.^ = e^,fi^) through a 6 ^ c 
quark transition, emitting primary leptons. 

3.1 Selection criteria 

A set of kinematic and quality criteria is used to separate hadronic events from 
continuum, QED and beam backgrounds, and to suppress contamination from cascade 
or fake leptons. 

Hadronic events are classified as having at least five tracks, a primary vertex with 
radial (K-) and z (Vz) components with respect to the origin satisfying V^ < 1.5 cm 
and \Vz\ < 3.5 cm, a total visible energy detected by ECL and CDC greater than 
50% of E^{^5) (^?(t5) being the energy of T(4S) at CM), a z component of the total 
visible energy less than 30% of -5^^(45) 5 ^ total calorimeter energy between 2.5% and 
90% of E'ym^n, and a ratio R2 of second and zeroth Fox- Wolfram moments [Q] less 
than 0.7. 

We select the electrons and the muons from the data sample of hadronic events. 
For the electron identification, we use the cluster energy, and position and shape of 
the electromagnetic shower in ECL, the dE/dx energy loss rate in CDC, the match- 
ing between the track in the CDC and the cluster in ECL, and TOF and ACC hit 
information for the fake lepton suppression. Optimization of the above parameters 
yields a ~ 90% efficiency for electrons and a ~ 0.3% misidentification probability for 
charged tracks with p > 1 GeV/c. Electrons consistent with pair production from 



photon conversions are removed. For the muon identification, we extrapolate the 
track in the CDC to the KLM, and we examine the difference between the expected 
and the actual penetration, as well as the matching quality of the KLM hits with 
the fitted track. The efficiency is ~ 85% for muon tracks and the misidentification 
probability ~ 2% for particles with p > 1 GeV/c. 

A "J/ip veto" is then applied, eliminating events consistent with a J/ip decay 
to a dilepton pair. All remaining leptons of the event are ordered according to the 
magnitude of their momentum at the T{4S) rest frame. Dilepton event candidates 
must have at least two particles identified as leptons. The leptons with the highest 
and the second highest momenta in the CM are selected as "dileptons" . Leptons are 
required to have a CM momentum and a laboratory polar angle within the ranges 
LI < p* < 2.3 GeV/c and 30° < ^'""^ < 135°, r and z distances of the decay 
vertex from the beam interaction point (BIP) satisfying |Ar^^^| < 0.05 cm and 
|Az^^^| < 2.00 cm, and at least two and one hits in the z and (r, 0) SVD layer strip 
planes, respectively. This is a combination of cuts aiming at suppressing cascade and 
fake leptons, as well as selecting events with good vertex resolution for the lepton 
tracks. To further reduce particle pairs originating from the same B, the opening 
angle between the lepton tracks calculated at CM, 9^^, is required to satisfy —0.80 < 
cos 9^^ < 0.95. The above values of the cuts have been optimized based on Monte 
Carlo (MC) studies. Application of all the selection criteria yields 8573 same-sign 
(SS) and 40981 opposite-sign (OS) dilepton events on the T(45') peak, and 40 SS and 
198 OS dilepton events below the resonance. 

3.2 Vertexing 

The proper time difference can be calculated from the z component of the distance 
between the decay vertices of the B meson pair, Az: At = Az/(3'~fc. The z- vertex po- 
sition of the leptons is determined by finding the intersection of each of the fitted lep- 
ton tracks with an ellipsoid known to contain the 5^ decay vertices. The dimensions 
of this volume are calculated by convolving the BIP "profile" (o"^^^ ~ 100 — 120 /im. 



ay ~ 5 /im, af ~ 2000 — 3000 fira) with the average flight length of the B meson 
(~ 20 /im at CM). The center and the dimensions of the BIP proflle are determined 
by the primary vertex position distribution on a run-by-run basis (a few tens of thou- 
sand events). The large uncertainty of the z component of the BIP is irrelevant since 
we are only interested in the difference between the z vertices, Az. Az is deflned as 
Zi+ — Zi- for the OS events, whereas for the experimentally indistinguishable leptons 
of SS events we use the absolute value of Az. 



Table 1: Classification of different terms in SS and OS dilepton samples and numbers 
of contributing events. 
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3.3 General strategy 

We classify the contributions to the SS and OS proper time distributions as "signal" 
(events with two primary leptons from semileptonic B^ or B^ decays), "background" 
{T{4S) events with at least one cascade or fake lepton) and continuum (all non- 
T(45') events). The background is further divided into two sub-categories indicating 
correctly or wrongly tagged events. The Az spectra for SS and OS dileptons can there- 
fore be described as sums of signal (S), correctly tagged background (C), wrongly 
tagged background (W) and continuum. Each of these contributions is characterized 
by a (normalized) proper time distribution, and is scaled by the number of source 
events (A^) and an overall selection efficiency (e). Parameters that control the popu- 
lations of source events are the number of T{4S) and continuum events in the data 
sample (A^t(4S) = N^^i ^^^ ^cont), the branching fractions of neutral and charged B 
pairs in T{AS) decays (/o and f± = ^ — fo), the fraction of mixed events in neutral 
B^ pairs {xd), and the semileptonic branching fractions for neutral and charged B 
mesons (6o and h±). Table |I| summarizes the categories of the events for the SS and 
OS distributions and the numbers of source events, A^. 

The signal terms for mixed {B^^B'^ / B^B'^: "mix"), unmixed {B'^B^: "unm") and 
charged {B^B^: "chd") B meson pairs are modeled analytically by smearing theo- 
retical expressions with explicit mixing and CPT-violating dependence: 



P„,ix(At) = i\sme\yATBo)e 



-|Atl/r. 



[l-cos(AmAt)] 



(1) 



Punm(At) = (l/4rBo)e '^*'^^^3[l + |cosep + (l- lcos^ncos(AmAt) 

-2/m(cos^)sin(AmAt)] (2) 

P,M(At) = (l/2rB±)e-l^*l/-i^± (3) 

A non-zero value for the complex parameter cos 9 would indicate CPT violation |Q . 
A cos 6' 7^ would also affect the B^ mixed event fraction, Xd '■ 

I sin 6*1^0;^ , , 

^'^ " |sin^|2x^ + (2 + x^ + x^|cos^|2) ^ > 

where Xd = r^o Anid ; 1/r^o = (Th + Tl)/2 is the mean width of the mass eigenstates, 

d d . . 

assumed to be known in this analysis |g]. 

If CPT is a good symmetry, then cos^ = 0, and the above expressions become 
the standard ones one finds in the bibliography: 

PLA^t) = (l/4r^o)e''^*'/^-^[l-cos(AmAt)] (5) 

PLJ^t) = (l/4r^o)e''^*'/^-^[l + cos(AmAt)] (6) 



Xd = 0/1 , ^2^ (7) 

^d) 



2(l + a;2^ 



The detector resolution for the signal terms is determined by the Az distribution 
of J/ip decays to lepton pairs, which here serves as a "response function". These 
leptons are selected with the same kinematic and quality criteria applied on primary 
lepton pairs from B decaysQ (Sec. fj.ll ). The Az distribution is obtained by finding 
two separate z vertices for the two lepton tracks of the J/ip decay, and by taking their 
relative difference. Its width is a = 112 /xm. To properly account for the long tails 
at large \Az\, the distribution is binned in a lookup table g, which is used to smear 
the theoretical functions: 

P(.. ^ ^ /^(At,^ear-At)P.(At)rf(At) 

"^ '"^'^^ //^(At,^ear-At)P„(At)rf(At)rf(Atsmear) ^^ 

with a = "mix", "unm", "chd". 

The modeling of the T(4S') background distributions is done numerically. To this 
end, we have simulated large MC event samples of generic charged and neutral B me- 
son Az distributions. The major component is the combination of a primary and a 
secondary lepton from a cascade decay (c — > s or r decay). We separate the contribu- 
tions from cascade and fake leptons. The number of cascade leptons is scaled to match 



-'^With the exception of the J/ip veto, of course. 



the B{B — > D"^ X), B{B — > D° X) branching ratios measured by CLEO 0. The 
number of fake leptons is independently determined with Ks — > tt+tt^ decays from 
the same run period. The small discrepancy in the detector's vertexing performance 
found between MC and data (~10%) is corrected by applying a convolution with a 
single Gaussian with a = 50ti2 /^m to each MC-determined background distribution. 
This procedure is tested on similarly obtained Az distributions from J/ip —>■ i^i~ 
and Ks -^ tt+tt" decays and is found to give a satisfactory matching. The strong 
dependence of the background distributions on the value of Arria is modeled through 
a linear interpolation of distributions from two generic B^ MC samples, generated 
with ArrLd = 0.423 ps~^ and Arrid = 0.464 ps^^. Finally, we have the option of ad- 
justing the steepness of the exponentially falling background distributions. This is 
useful when we want to perform the fit for a different value of the (fixed) B^, B^ 
lifetimes, or calculate the associated systematic error. 

The continuum Az distributions are also modeled numerically with MC simulation 
samples. The relative weight of the continuum contribution is determined by the 
number of off-resonance data, scaled to account for luminosity and energy differences. 
The correction to the vertex resolution is applied here as well. 

3.4 Fitting 

The data is fit to a linear combination of analytical and numerical expressions to 
extract the mixing and the CPT violation parameters. A binned maximum likelihood 
fit is performed simultaneously to the SS and OS dilepton Az distributions in the 
region \Az\ < 1.85 mm Q]. 

We fix the parameters r^o = 1.548 ps §, TBi/r^o = 1.06 §, /±//o = 1.05 § 
and 6±/6o = 'Tb±/'Tb° 01 • Besides Arrtd and the (complex) parameter cos 6', the fit 
includes two more free parameters, the selection efficiency ratios e^/ec and e^/ec for 
the charged B mesons. The relative selection efficiencies for mixed, unmixed and 
charged B mesons within each event tag type (S*, C or W) are fixed to their MC 
values. 

The expected total MC distributions for the SS, OS spectra are not normalized 
separately to the numbers of the events of the two histograms, but only to the total 
number of data events Ntot = Nss + Nqs- In other words, the ratio of the populations 
of the SS and OS total MC distributions is determined by the fit result for the mixing. 

Assuming that CPT is a good symmetry, and therefore by fixing cos^ = 0, we 
find Arud = 0.463 ± 0.008 ps"^ with x^/Ndoi = 332.5/376. The efficiency ratios are 
es/ec = 143.2 ± 5.1 and es/ew = 103.6 ± 15.0. 

When cos^ is a free parameter in the fit, we find Arrid = 0.461 ± 0.008 ps~^, 
Re{cose) = 0.00 ± 0.15 and /m(cos^) = 0.035 ± 0.029, with x^/Ndd = 331.0/376. 
These results are consistent with CPT symmetry. The efficiency ratios in this case 



are es/ec = 142.8 ±5.1 and es/^w = 102.8 ± 14.9. The fractions of signal events 
are found to be 32.1% in SS (mixed B pairs) and 77.5% in OS (32.1% unmixed and 
45.9% charged B pairs). 

Fig. ID shows the A2; distributions for the data and the fitted Monte Carlo distri- 
butions, for the SS (top) and OS (bottom) dilepton spectra. The signal in the SS 
histogram is a small term compared to the background from B^B^^ and B^B^ . A 
much more evocative way of displaying the B^ — 5° oscillation is by plotting the time 
dependent asymmetry between SS and OS dileptons: 

,,, , Ar+-(At)-iV±±(At) ,^ , , 

M^t) = .r^ /A .r^^)^ ~ COS Am At 

^ ^ iV+-(At) + A^±±(At) ^ ^ 

To a good approximation, most of the background cancels out, and the time 
dependent asymmetry reveals a clean oscillation signal, as illustrated in Fig. 0. In 
the same plot, a convolved cosine term is superimposed on the data points, with the 
period determined from the fit. 

3.5 Systematic errors 

To estimate the systematic errors of the various assumptions made in the fit, the 
associated input parameter is varied, typically by ±a, and then the fit is repeated. 
The shift in the central value is taken to be the systematic error for that category. 

The main contributions to the systematic errors are from the uncertainty in the 
measurements of the B lifetimes and the branching fractions of T(4S') to charged 
and neutral B pairs. We calculate these contributions by adjusting the corresponding 
parameters by the amount of their uncertainties, one at a time. The lifetime changes 
affect all 5-originating terms in the fit (exponentially decaying distributions) and, in- 
directly, signal or background events with one primary lepton through the propagated 
change in the semileptonic branching ratios. The change of branching fractions f±, 
/o affects the relative populations of charged and neutral B pairs in the distributions 
(signal and background). 

We have made the approximation that the difference between the decay widths 
of the 5° mass eigenstates, AF^, is zero. The conservative estimation used by most 
phenomenologists is AF^/Frf < 10~^ 0. We report the shift in the fitted values by 
assuming a ATd/Td = 1% and by modifying accordingly the proper time distribu- 
tions 0]. For the systematic errors associated with the response function, we estimate 
the effect of the B motion in the T(4S') frame (which includes the discrepancy of the 
J/ip A^ distribution from the true dilepton response function) and the error originat- 
ing from the statistical uncertainties in the determination of the response function. 
This set of changes only affects the signal terms. 
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Figure 1: Az distributions for same-sign (top) and opposite-sign (bottom) dileptons 
for data (points) and the fitted Monte Carlo distributions (histograms). For the 
Monte Carlo distributions, signal and background from different categories are plotted 
separately. 
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Figure 2: Time dependent asymmetry between SS and OS dileptons. The curve is a 
convoluted cosine with period determined by the fit. The negative A2; region for the 
OS histogram has been folded into the positive region for display purposes. 
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We change the cascade lepton rates by modifing the branching ratios of B mesons 
decaying to D^ and D^ in the MC by ±4.6% and ±14.3%, respectively. We change 
the levels of fake lepton and continuum contamination in the sample by shifting the 
fake rates by ±35%, and the continuum contribution by ±16 (SS) and ±7% (OS). The 
error in the determination of the width of the convolving Gaussian (a = SOlJg fim) is 
used as the uncertainty in the detector's vertexing performance for the background 
terms (~ 30% of the sample). We change the width of the Gaussian and report 
the shift in the fit central values with the new background distributions. We check 
the effect of finite MC statistics by varying the number of background events on a 
bin-by-bin basis by the amount of the statistical errors. 

We finally check the effect of binning and several other possible sources of system- 
atic errors, and we find that their contribution is negligible. 

Table 1^ summarizes the systematic errors for Am, Re{cos6) and Im{cos6). Since 
the asymmetry between positive and negative errors is small (~ 15 — 20%), we con- 
servatively take the larger of the values to be the systematic error. 

4 Summary 

We have measured the mixing parameter Arrid and searched for CPT violation in 
the B^ system from the time evolution of dileptons in T(45') decays. The analysis 
corresponds to an integrated luminosity of 5.9 fb~^ of data collected with the Belle 
detector. The run period was from January to July of 2000. 
If we invoke CPT symmetry for Am^ we obtain: 

Arrid = 0.463 ± 0.008 (stat) ± 0.016 (sys) ps"^ 

If we fit simultaneously for Arria and the CPT- violating parameter cos 6 we obtain: 

Arrid = 0.461 ± 0.008 (stat) ± 0.016 (sys) ps"^ 
Re{cos9) = 0.00 ±0.15 (stat) ±0.06 (sys) 
Jm(cos^) = 0.035 ± 0.029 (stat) ±0.051 (sys) 

These results imply [Q for the mass and the lifetime differences between i?° and B^ : 

Im^o -m^ol/niBO < 1.6 x 10"^^ and 11^0-1^01/1^0 < 0.161, at 90% C.L. 

d d d d d d 

The results are consistent with CPT conservation. This is the first direct measure- 
ment and the first published result of Amd from a time-dependent analysis on the 
T(4S') resonance p|. The limit on (F^o — r^o)/r^o measurement is consistent with 

d d d 

previous measurements p. This is the first measurement on Re{cos6) and on the 
mass difference between B^ and -B°. 
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Table 2: Summary of systematic errors (a) with and (b) without invoking the CPT 
symmetry. 
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Fractions of charged and neutral B meson pairs 
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+0.007 
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Continuum 
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Detector resolution {tH /im) 


+0.001 


-0.001 
+0.000 
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-0.000 


MC statistics 
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±0.005 
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±0.009 


Other systematic errors 


Binning 


±0.001 


±0.001 


<0.01 


±0.002 


Total 


+0.012 
-0.016 


+0.014 
-0.016 


±0.06 


±0.051 
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